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Abstract 

A solid sphere is considered, with a uniformly distributed infinity of points. 
Two points being pseudorandomly chosen, the analytical probability density that 
their separation have a given value is computed, for three types of the underlying 
geometry: E^, and S^. Figures, graphs and histograms to complement this short 
note are given. 
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1 Introduction 



The shape of the universe is presently under investigation, and cosmic crystallography 
(CC) is one of the various methods proposed to determine it [jll. 

The idea which supports the CC method is that if the universe is multiply connected 
then multiple images of a same cosmic object (a given quasar, say) may be observed in 
the sky. The separations between pairs of these images are correlated by the geometry 
and the topology of the spacetime; so if one selects a catalogue of observed images of 
cosmic objects and performs a histogram of the separations I between the images, then 
these correlations manifest either as spikes (associated with Clifford translations) or as 
slight deformations of the histogram of the corresponding simply connected manifold 0]. 

A variant method was recently proposed by Fagundes and Gausmann , of subtract- 
ing from a histogram 0(/j) of a multiply connected space a histogram of the corresponding 
simply connected space. The reported result was a plot with much oscillations in small 
scales. 

In the present paper we propose an alternative to Fagundes- Gausmann method: we 
derive continuous probability density functions to be subtracted from the histogram 
0(/j) , and thereby obtain a histogram sensibly more suitable for analysis. 

In section 2 we derive the functions jF(a, /) for the Euclidean, hyperbolic and elliptic 
geometries, and in section 3 we make a few comments. 

2 Probability densities 

We consider one of the simply connected spaces E^, or S^. In that space, a spherical 
solid ball Ba is taken , with radius a. The ball is assumed to contain an infinite number 
of pointlike objects, spatially distributed as uniformly as possible. We next select pseu- 
dorandomly two points of Ba and ask for the probability jF(a, l)dl that the separation 
between the points lie between / and I + dl. 

a. Euclidean geometry 

In the Euclidean three-space take a ball Ba centred at the origin O and select two 
points P, Q in the ball; let r G [0, a] be the radial position of P and let / < 2a be the 
distance from P to Q (see Figure 1). 

Clearly the probability density J-'e{0', r, I) of this configuration is proportional both to 
the area Ssir) = 47rr^ of the locus of P well as to the area of the locus of Q. When 
r + I < a this locus is a sphere S"^ with area Se{1), while when r + I > a the locus is a 
spherical disk in with area DE{ci,r,l) = (7r//r)[a^ — (Z — t)^]. 
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The probabihty density of the configuration is then 



J^E{a,r,l) = kSE{r){SE{l) x Q{a - I - r) + DE{a,r,l) x Q{l + r-a)}, (2.1) 

where k is a normahzation constant and G is the Heaviside function. 
We integrate eq. ( ^.11 ) for r G [0, a], and finaUy obtain, for I G (0, 2a], 

^E{a,l) = ^{2a-l)\l + 4a), (2.2) 
where the value k = 9/{167['^a^) was set to satisfy the normahzation condition 
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J^E{a,l)dl = l. (2.3) 



One often encounters in the hterature reference to the probabihty density V{s) that 
the squared separation be s; since s = f and V{s)ds = J^{l)dl, then 



VE{a. s) = |^(2a - ^sf{^s + 4a), (2.4) 



vahd for s G {Q,Aa?]. 

In Figure 2 we reproduce a typical mean pair separation histogram (MPSH) for pseu- 
dorandomly distributed objects in an Euclidean solid sphere Ba with arbitrary radius, 
together with the corresponding probability density ^^^(a, /). It should be stressed that, 
differently from the hyperbolic and elliptical cases, the shape of the function TEifl-, does 
not depend on the value of the radius a. 

b. Hyperbolic geometry 

To obtain the probability density Tuip^^ ^^e hyperbolic geometry we follow the 
same lines as before. The area of a sphere with radius r is now Siiir) = 47ri?2 sinhV/i?, 
where R is the radius of curvature of the geometry; without loss of generality we henceforth 
set i? = 1. On the other hand, the area of a spherical disk in is (see Figure 1) 

Dnia, r, I) = 2tt sinhZ[ sinh/ — cosh/ cothr + cosha cschr], (2-5) 

to be considered whenever r + / > a. 

The probability density of the configuration is then 

jr^(a,r,/) = kSH{r)[SH{l) x e{a - I - r) + Dnia^rJ) x e(/ + r-a)], (2.6) 

which upon integration for r G [0, a] and normalization gives 

8 sinh^ / 

jr^(a,/) = - — [cosha sech(//2) sinh(a - 1/2) - (a - 1/2)], (2.7) 

(sinh2a — 2a)^ 
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valid for Z e (0,2a]. 

In Figure 3 we reproduce graphs of J^nia, I) for three values of the radius a. For a « 1 
the function tends to the Euchdean one given in Figure 2, as expected. For increasing 
values of a the function shifts towards the large values of I, and for a » 1 a strong 
concentration of Th{ci, I) is found near the extreme value I — 2a. 

Figure 4 shows a typical MPSH in the hyperbolic three-space, together with the cor- 
responding probability density Th{(i, 0- 

c. Elliptic geometry 

The basic strategy to obtain the probability density J-'s{a,l) is the same as before, 
and the calculations are similar whenever the diameter 2a of the ball is less than the 
separation nR between antipodal points in the three-sphere; we then find for J-s{ci, I) the 
expression, valid for / e (0,2a], 

8 sin^ I 

Ts{a < 7r/2, /) = — . [(a -1/2)- cosa sec(//2) sin(a - 1/2)], (2.8) 

yZiCL SlH ^CL) 

where we have taken R—1 without loss of generality. 

However, the cases where a > 7r/2 are considerably trickier to deal with, due to the 
treacherous connectivity of the spherical space and the new requirement that I must 
not exceed tt. A much larger quantity of trivial integrations now comes into scene, before 
the following expression is eventually obtained: 

8 sin^ I 

= |2a-sm2ap »^°-^'°^''-^'+ f^-") 

©(27r — 2a — /) x [sin2a -\- -k — a — 1/2 — cosa sec(//2) sin(a — 1/2)]}, 

vahd for all a e (0, tt] and I e (0, min(2a, tt)]. 

In Figure 5 four graphs of !Fs{a, I) are shown, for different values of the radius a of 
the ball. For a increasing from to tt the function continuously shifts towards the smaller 
values of L In particular, when a = 7r/2 we have 

Ts{t^/2, /) = -(!- //tt) sinH , (2.10) 

TT 

while when a = tt we have the harmonic, symmetric probability density 

J^s{'K,l)^-smn . (2.11) 

TT 

In Figure 6 a sample MPSH in the spherical space is reproduced, together with the 
corresponding probability density ^5 (a, /). 
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3 Comments 



It is perhaps worth clarifying the meaning of the functions jF(a, r, /): if we pseudorandomly 
choose two points P and Q in a sohd sphere with radius a (see Figure 1) then jF(a, r, l)drdl 
is the probabihty that P hes between the radial positions r and r+dr, times the probability 
that the separation from Q to P lies between the values / and l + dl. The form of jF(a, r, /) 
clearly depends on the geometry one is concerned with. 

Each histogram in Figures 2, 4, and 6 has m = 100 subintervals and is a MPSH - a 
mean pair separation histogram comprising K = 10 comparable catalogues with = 50 
objects each All computer-generated histograms assume a homogeneous distribution, 
as described by Lehoucq, Luminet and Uzan [Q. 

A close inspection of Figure 2 shows that the most probable separation between two 
arbitrarily chosen points in an Euclidean solid ball is slightly greater than the radius of the 
ball; also the maximum of Ve{0'iS) in eq.( |2.4D occurs when s/4a^ = 0.134, in agreement 
with plots of Fagundes and Gausmann p. 

A characteristic feature of the hyperbolic geometries is that at large distances there 
is more space than in the Euclidean geometries; this fact is clearly exhibited in Figure 3, 
which shows a strong predominance of large separations I when the radius a of the solid 
sphere is large. It is worth noting that Fagundes and Gausmann |^ obtained histograms 
with a = 2.34 which closely resemble ours with a = 2 in Figure 4. In contrast, the 
hyperbolic histograms given by Lehoucq et al. ^ use radius a nearly 0.6, so they are 
similar to our Euclidean one. 

Oppositely to the hyperbolic case, in distant places endowed with the elliptic geometry 
there is less space than in the Euclidean geometry; this is nicely illustrated in Figure 5. 
Indeed, with increasing a ( increasing solid ball) the probability of finding small distances 
/ (in comparison with a) also increases. We further note that when a increases from 7r/2 
to 7r the ratio Imax/'^O' recedes from 1 to 1/2. 
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Captions for the Figures 



Figure 1. A two-dimensional picture of the solid ball Ba- We have OR = a (the radius 
of Ba), OP = r, PQ = I. The circular arc with centre P represents a spherical disk D^; 
the disk becomes a sphere S"^ whenever r + I < a. 

Figure 2. In the Euclidean space E^, a sample MPSH for a solid ball with arbitrary 
radius a, together with the corresponding probability density eq. (|2.2|). 



Figure 3. In the hyperbolic space with curvature radius R = 1, the probability 
densities eq.(|2.7]) for solid spheres with radii a = 0.01 (dots), 2.0 (line) and 10 (bold line). 



Figure 4. A sample MPSH and the corresponding probability density eq. (|2.7|) for a 
solid sphere with radius a = 2 in with unitary radius. 

Figure 5. In the three-sphere with radius R = 1, the probability densities eqs.( |2.9| )- 
( p.ll ) for balls with radii a = vr (bold line), 37r/4, 7r/2, and 0.01 (dots). 



Figure 6. In the spherical space with unitary radius, a sample MPSH for a solid 
ball with radius a = n/2, and the corresponding probability density eq. (|2.10|) . 
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